 use "$data\women_analysis.dta", clear
  
 global controlwoman mbase_ageyears  fdstotal rantime ranerror numberkids
 global controlhh land members adults mbase_electrified  mbase_OB mbase_scheduled
 

 estimates clear
 foreach i in confidencescoregeneral mobility_index_do Pleave1 Pmarketin1 Pmarketout1 ///
 Phealth1 Pnatal1 Pcallnatal1 Pgroup1 Pmela1{
 ivregress 2sls `i' (participatedTA_IV=ITTanalysis) i.strata $controlwoman $controlhh if sequence==1, vce(robust)
 est store `i'
 su `i' if sequence==1 & ITTanalysis==0
 estadd scalar controlmean=`r(mean)'
 }
 
 esttab * using "$results\Table8.xls", stat(controlmean N r2) keep (participatedTA_IV) b(3) se(3) replace
 
 ***FDR q-values
 *The FDR Q-values are computed following  Benjamini and Hochberg (1995) within each category.
 
 exit

